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1. Introduction 

The search for the Higgs has become a major issue in particle physics as the LHC is nearing 
its completion. The Standard Model (SM) cannot be considered complete given that the Higgs 
is as-yet unobserved and it is not clear how Electroweak symmetry is broken in nature. If the 
Higgs is seen, its properties could tell us about physics beyond the Standard Model, such as the 
energy scale of a more fundamental theory. The current lower bound for the Higgs mass from 
direct searches is 114.4 GeV [|l]]. The Higgs mass can also be inferred indirectly by fitting the 
Standard Model to a host of Electroweak precision measurements. The best perturbative fit gives 
niH = 76^24 GeV, so the data certainly seem to prefer the Higgs to be light However, the global 
fitting procedure, which favors a surprisingly low Higgs mass, has its own intrinsic issues, perhaps a 
hint that deviations from the Standard Model are already present [^]. Larger Higgs masses together 
with new physics threshold effects at the TeV scale will require new extended analysis |Q, ||] where 
non-perturbative effects may come into play. 

Based on the assumption that the Standard Model is only valid up to some energy scale, lower 
and upper bounds on the Higgs mass were established before without relying on input from Elec- 
troweak precision measurements. Bounds on the Higgs mass are valuable for two reasons. Firstly, 
they cut down the parameter space where one searches for a Standard Model Higgs. Secondly, if 
the Higgs is found, measuring its mass and knowing the bounds it must obey would indicate the 
maximum energy scale up to which the Standard Model can work. In phenomenology, the origin 
of the lower bound is thought to be the vacuum instability the Top quark loop would generate, if 
the Higgs mass were too light. The upper bound in phenomenological analysis is simply calculated 
by not allowing the running Higgs coupling X{t) to become strong at the cutoff scale A which 
represents new physics before X{t) would run into the fictitious Landau pole. These ideas on lower 
and upper Higgs mass bounds have been applied to the Standard Model for almost 30 years and 
have been increasingly refined. 

The bounds given by the state-of-the-art cal- 
culations were reviewed in [^] and shown in Fig- 
ure H, based on the original work in [^] and [||]. 
There are several things one can learn from this 
plot. The Standard Model apparently cannot gen- 
erate a Higgs boson heavier than 1 TeV without 
strong Higgs self-interactions and a low thresh- 
old for new physics in the TeV range, a scenario 
not consistent with the perturbative loop expan- 
sion of the Electroweak precision analysis. What 
non-perturbative modifications on the TeV scale 
would support a heavy Higgs particle, consistent 
with Electroweak precision data, is one of the mo- 
tivations for our lattice studies [Q, |5|]. The lower 
bound is interesting for today's phenomenology, 
given the current experimental limits. If the Higgs 
mass is around 100 GeV, this would intersect with the lower bound in Figure |l| somewhere between 
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Figure 1: Upper and lower bounds for the Higgs 
mass as a function of the scale of new physics be- 
yond the Standard Model, from [|]). 
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10 and 100 TeV, beyond which apparently new physics should enter. 

One major goal of our lattice Higgs project is to understand the role of vacuum instability 
and the Landau pole in an exact non-perturbative setting when the intrinsic cutoff in the Higgs 
sector is not removable and low in the TeV range. Another goal is to explore the role of non- 
perturbative Higgs physics from the lattice in extensions of the perturbative SM analysis, including 
the possibility of a heavy Higgs particle within the Higgs reach of the LHC. 

The outline of this paper is as follows. In section 2 we will report results from the large 
Np analysis of the Top-Higgs Yukawa model of a single real scalar field coupled to Nf fermions. 
The influence of the non-removable intrinsic cutoff (triviality) on the exact renormalization group 
(RG) flow is exhibited. The vacuum instability problem of the model is discussed on the lattice 
in section 3 and compared with the traditional renormalization group procedure of the Standard 



Model (earlier versions of this work on vacuum instabihty have been discussed in [g] and []10|]). 
In section 4 we present the Wilsonian view on the renormalization group as applied to the vacuum 
instability and Higgs lower bound problems. The first lattice simulation results on the Higgs mass 
lower bound, using chiral lattice fermions in Top-Higgs Yukawa models, are reported in section 5. 



Using the higher derivative (Lee-Wick) extension of the Higgs sector Jl 1| , |12| , |13[ ], we will 
illustrate in section 6 how non-perturbative lattice studies might help to investigate heavy Higgs 
particle scenarios in the 500-800 GeV Higgs mass range relevant for future LHC physics. Con- 
straints from Electroweak precision data on the heavy Higgs particle are briefly discussed. 

2. Top-Higgs Yukawa model in large Np limit 

For pedagogical purposes, we first consider a Higgs- Yukawa model of a single real scalar field 
coupled to Np massless fermions. The saddle point approximation in the large Np limit becomes 
exact and this will allow us to demonstrate that the theory is trivial. We will also calculate the 
flow of the renormalized couplings as a function of the energy scale to identify problems with 
the vacuum instability scenario when the intrinsic cutoff is non-removable. Similar behavior is 
expected at finite Np which requires non-perturbative lattice simulations. 

2.1 Renormalization scheme 

Let us start with the bare Lagrangian of the Higgs- Yukawa theory in Euclidean space-time, 
which is 

^ = + ^Ao0o' + \ {d^hY + wS {r^d^ +yo(h) vS, (2-i) 

where a = l,...,Np sums over the degenerate fermion flavors and the subscript denotes bare 
quantities. We rewrite this as 

+ (1 + 5z^W7nd^,V + + (2.2) 

where we have introduced the wavefunction renormalization factors = 1 + 5z0 , Z^, = 1 + 5z\ir 
and renormalized parameters with their corresponding counterterms. The connections between the 
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bare and renormalized parameters are 

mpZ^ =n? + 5m^, AqZ^ = A + 5A , Z^y^yo =y + dy. (2.3) 

In the limit where Nf becomes large, all Feynman diagrams with Higgs loops are suppressed rel- 
ative to those with fermion loops. Hence, two of the counterterms vanish, 5y = 0, Szxf, = 0, as 
there are no radiative corrections to the fermion propagator or to the Higgs-fermion coupling. Let 
us specify the renormalization conditions which determine the remaining counterterms. 
In the large A^^ limit, the renormalized Coleman-Weinberg effective potential 1 14 ] is 

C/eff = ^m202 + _LA0V^5m2(^)2 + i-5A04-2Aff /ln[l+/0V^'] (2.4) 

Z Z4 Z Z4 Jk 

containing the tree-level contributions from the renormalized parameters and their counterterms, 
and the infinite sum of all diagrams with one fermion loop and an even number of external legs. 
The factor Np comes from all the possible fermions which can appear in the single loop and we use 
the notation J d'^k — > Jj^ for loop integrals. The vacuum expectation value (j) = v is where Usff 
has an absolute minimum i.e. ?/gg(v) = 0. In the Higgs phase of the theory, v 7^ 0. At tree-level, 
this gives the relation 

+ \Xv^ = 0, (2.5) 



coming from the first two terms in Equation (2.4). Our first renormalization condition is that we 
want to maintain the tree-level relation in Equation ( ^ ) exactly, giving 

5m^ + ^dXv^ - ANpy^ I \ , = 0. (2.6) 

The counterterms exactly cancel all the finite and infinite contributions of the radiative diagrams. 
The same relation can also be determined by demanding that the tadpole diagram is exactly can- 
celled by the counterterms. 

In the Higgs phase, we define the Higgs fluctuation around the vev as = (p + v. At tree-level, 
the mass of the Higgs fluctuation i.e. U'^^{v) is 

= + ^Av^ = ^Av^ (2.7) 
In the large A^^ limit, the inverse propagator of the Higgs fluctuation is 

G-^lip") = p^ + m^+^-Xv^+p^5z^ + dm^ + ^5Xv^-Up^) 

= -W,y' I ...^"^w^'^^/^l 2 2V (2.8) 
Jk [k^ +y'^v'^)[{k- pY +y^v^) 

where all Higgs-loop diagrams are suppressed relative to the single fermion-loop diagram. We 
impose the condition that 

G-^lip" ^0)=p'' + ml, (2.9) 
which separates into two renormalization conditions: 

5m^ + ]^8Xv^-I.{p^ = 0)=0 (2.10) 
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and 



5z0 



0. 



(2.11) 



p^-=0 



The renormalization condition Equation ( |2.10 ) maintains the tree-level relation in Equation (2.7) 
exactly. Again, the counterterms precisely cancel all the finite and infinite radiative contributions. 
We should point out that the Higgs mass defined as the zero-momentum piece of G^^ is identical 
to that defined via the curvature U^ff{v). This is not the same as the true physical mass given by the 
pole of the propagator, and these masses can be related to one another in perturbation theory. 

The renormalization conditions Equations (2.6) and ( 2.10| ) can easily be solved. Because we 
wish to demonstrate triviality in this theory, we use some finite cutoff in the momentum integrals 
and examine what occurs as this cutoff is removed. We will use a simple hard-momentum cutoff 
1^1 < A. Exactly the same conclusions would be reached using instead e.g. Pauli-Villars regular- 
ization. The non-zero counterterms after the loop integration are 



6m^ 
5z,j> 



2712 

3NFy 



,4 



4 4 

y V 



1 



2 2 



2 2 

y V 



1 1 



In 



2 2 

y V 



2(A2+3;2v2) 2 2 ^y^v'^ 
3;V+A2^^ ^ -5a4 _ 3A23;2v2 
12(A2+3;2v2)2 



Iln 



3;2v2 



(2.12) 



As we said earlier, in the large Np limit, the fermion inverse propagator receives no radiative 
correction, 

G^\[P) = P\iy\^ +3'v, (2.13) 

so we identify the fermion mass as mt = yv (looking ahead to the Top quark), which we substitute 
into all of the above equations. 

2.2 Triviality 

Let us first consider the regime mj/A ^ 1, where the cutoff is much larger than the physical 
scale. In this limit, we get 



1 + 



8712 



In 



rA2i 




2 

rUf 





(2.14) 



For any finite bare Yukawa coupling yo, the Higgs wavefunction renormaUzation factor vanishes 
logarithmically as the cutoff is removed, nit/K^O. This same logarithmic behavior, for any choice 
of bare couplings, will appear in all of the renormalized couplings, leading to the triviality scenario: 
a finite cutoff must be kept to maintain non-zero interactions. Explicitly, the renormalized Yukawa 
coupling is 



y^ = ylZtj, = yl 



1 + 



8712 



In 



A 



21 



mf 



-1 



^In — 
87r2 m} 



as 



(2.15) 
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For the renormalized Higgs coupling, we have 



AflZ^ — 5X 



,4 / 1 



1 mf\ 

2'" 7?) 



2(A2+m2) 
2- 



2 2"VA2 + m2 



12 



-1 



as 



nit 



0. (2.16) 



The slow logarithmic vanishing of y and A allows to have a relatively large separation of cutoff and 
physical scales and still maintain significant interactions. However, the standard renormalization 
procedure of removing the cutoff completely gives a non-interacting theory. Although completely 
unphysical, we can also consider the limit m, /A » 1, where the cutoff is much below the physical 
scale. From Equation (2.12), we see this gives 5X = 0, 5z0 = 0, and hence Z^ 1. In this limit, 
the connection between bare and renormalized parameters is simply A = Aq, j = jo- This result is 
not surprising: deep in the cutoff regime, we simply have the bare theory, with no separation into 
renormalized parameters and their counterterms. This will be relevant when we discuss whether 
the vacuum can become unstable. 



2.3 Renormalization group flow 

The physical properties of the theory are fixed as soon as one chooses a complete set of bare 
parameters. As the cutoff is varied, the renormalized couplings flow in order to maintain exactly 
the renormalization conditions we have imposed. Using the explicit cutoff dependence of y and 
A, we can calculate this Callan-Symanzik flow. In the limit m(/A <C 1, from Equations ( 2.15D and 



(2.16), we have 



A 



A 



1 



' 4712 



dA 16712 



[-8A^fA/+48A^f/ 



(2.17) 



The same /3 functions would be obtained in the large Np limit for the running y and A couplings 
in scale dependent RG flows using e.g. dimensional regularization, where no cutoff would ex- 
plicitly appear. (Since increasing A corresponds to decreasing mass scale jj., the j3 functions in 



Equation (2.17) have opposite signs). It is important to note that the two RG schemes have very 



different physical meanings: Equation (2.17) describes the response to changing the cutoff whereas 
the scale dependent RG flow compensates for the arbitrary choice of the renormalization scale at 
finite cutoff. When the cutoff is far above the physical scales, the finite cutoff effects are neg- 
hgible and we expect to reproduce the unique cutoff-independent j8 functions. However, as the 
cutoff is reduced and m,/A increases, this cannot continue to hold indefinitely, as the renormalized 
couplings must eventually flow to the bare ones, as explained above. 

Let us demonstrate an explicit example of the Callan-Symanzik RG flow in the presence of a 
finite cutoff. In the large Nf limit, nit =yv = yoVQ. The bare vev is determined by the minimum of 
the bare effective potential 



1 2 2 1 



-2Nf y^ln [1+3^202/^2 



(2.18) 
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0- 



m^p/cutoff m^/cutoff 



Figure 2: The exact RG flow of the renormalized couphngs X and y with the fuU cutoff dependence. The 
corresponding bare couplings are Aq = 0.1 and yo = 0.7. For large cutoff, the exact flow agrees with the 
continuum RG flow, where the cutoff dependence is omitted. For small cutoff, the exact RG flows to the 
bare couplings Aq and yo, but the continuum RG misleadingly predicts that X turns negative. 



Using a hard-momentum cutoff, this gives 



ml + ^hvl ■ 



27r2 



(2.19) 



We express all dimensionful quantities in units of the cutoff A. We pick some fixed values for 



Ao and 3^0- Varying the value of htq/A^ changes the solution vq/A of Equation ( 2.19 ) and hence 
the ratio m,/A. As we said, choosing the values of the bare parameters completely determines 
everything in the theory. For example, to attain a very small value of nit /A requires nig/A^ to be 
tuned quite precisely. Using Equation ( 2.19| ), the critical surface, where vq/A = 0, is the transition 
line 



A2 



47r2 



0. 



(2.20) 



Using Equations (2.3) and (2.12), all of the counterterms and renormalized parameters can be 
expressed in terms of ?io, yo,ml and vq- Solving this set of simultaneous equations is a simple 
numerical exercise. We make an arbitrary choice Ao = 0.1, = 0-7 which would correspond to 
the physical Higgs below its lower bound in phenomenological considerations. Varying the value 
of ml/A^, we explore numerically the range 10^'^ < m^/A < 102. The results in a limited range 
are plotted in Figure 0. When the cutoff is high, the exact RG flow is exactly the same as if the 



cutoff had been completely removed and follows precisely the continuum form of Equation ( |2. 17| ). 
However, as the cutoff is reduced, the exact RG flow eventually breaks away from the continuum 
form and reaches a plateau at the value of the bare coupling. 

The continuum RG in the above example predicts that A turns negative at some energy scale 
as the flow continues. This was used in the past as an indication that the ground state of the theory 
turns unstable at that scale which would determine the energy scale of new physics necessary to 
sustain a particular value of the physical Higgs mass (vacuum instability bound). As shown above, 
the true RG flow with the full cutoff dependence saturates at Ao and does not turn negative under 
the necessary Ao > stability requirement of the model. This makes the phenomenological RG 
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method and the apparent vacuum instabiUty quite suspect in the presence of the non-removable 
finite cutoff which is required by triviality of the renormalized couplings. 

The absence of vacuum instability will be demonstrated directly in the next section using the 
Higgs effective potential. In sections 4 and 5 we will propose a lattice strategy to determine the 
Higgs mass lower bound in the presence of an intrinsic cutoff without relying on the continuum 
RG flow. In this new strategy even the Ao > condition might be relaxed by adding new irrelevant 
operators, like the ^<^^ term, to keep the stability of the cutoff theory intact. 

3. The effective potential and vacuum instability 

First, we will present here the RG improved one-loop calculation of the effective potential 
with unstable vacuum when the cutoff is ignored. Next we show the absence of vacuum instability 
when the cutoff is correctly enforced. 

3.1 Continuum 1-loop effective potential 

For the Higgs-Yukawa model with Np fermions of Section 2, the 1-loop renormalized effective 
potential is 

f/efif = ^mV + ^A0V^W02 + ^5A04-2A^f^ln[l+jV/^2] 

V = \ni^^^ + l-x^\ (3.1) 
where the Higgs-loop contributions are also included now. For consistency, we impose exactly 



the same renormalization conditions Equations (2.10) and (2.11) used in Section 2, including all 



the Higgs-loop radiative corrections. Because 5y and 5z\^ are non-zero (we no longer impose the 
large Np limit), we specify the two additional renormalization conditions. The fermion inverse 
propagator is 



G^^\,{p) =p^ly^+yv + 5zx^Pn 7^ + ^3^^ - If (/?) , 
2 f -k^Yn+yv 

h{B+y2v2){{k-pf+ml 



the radiative correction coming from a single Higgs-loop diagram, and we require that 

G^,\f{p^^)=Ptiy^+yv. (3.3) 

This gives two renormalization conditions, 

5jv-rf(/7^o) = , 



5z 



(3.4) 



Again, the counterterms completely remove all the finite and infinite parts of the radiative correc- 
tions. We regulate the momentum integrals using e.g. a hard-momentum cutoff. The counterterms 
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and the renormalized effective potential are calculated exactly using a finite cutoff. We then take 
the naive limit 0/A ^ to remove all cutoff dependence. This ignores the fact that a finite and 
possibly low cutoff is required to maintain A,j / (a crucial point why the instability does not 
occur in the presence of finite cutoff). 

The continuum form of the 1-loop renormalized effective potential is given by 



2 ^ 24 ^ 16712 
1 



-0V2v^0" + 0^1n^ 
2 



167r2 

+ 



1 n2w.4 .2 2^, ^2 + ^02/2 

— {^ r - l^rmji) In ^—^ 

lo mi 



H 

1 4 m2 + A02/2 3 ,2a4 , "7 1 a2 2 

16 32 16 



(3.5) 



where = Av2/3. Due to our choice of renormalization conditions, the tree-level vev v 



3mfj/X is not shifted: one can check explicitly that U^s in Equation ( p^ ) has its minimum 



at = V. The large Np limit can be recovered by omitting the Higgs-loop terms. 

3.2 RG improved effective potential and vacuum instability 

The stability of the ground state is determined by the behavior of Usn for large (j). We see 
from Equation (3^) that the dominant terms in this regime are of the form A^0'^ln(02/v2) and 



—Nfy^(j)'^ln{(l)^/v'^). The negative fermion term brings up the possibility that the vev v is unstable. 
Hence stability is determined by the relative values of A 2 and y^, which are related to m// and m,. If 
the fermionic term dominates at large (p, the minimum at v is only a local one and will decay. If we 
believe that the vacuum is absolutely stable, then new degrees of freedom must enter at the scale 
where ?7eff(0) first becomes unstable. For given values of mn and nit, this predicts the emergence of 
new physics. Turning this around, let us fix m, and ask that no new stabilizing degrees of freedom 
are needed for < Then we obtain a lower bound mu{E): if the Higgs is lighter than this, U^s 
is already unstable for below E because the fermion term dominates even earlier. 

Improved vacuum instability can be shown via the running renormalized couplings in RG 
setting. We can define a set of renormalization conditions in the continuum, for example in the MS 
scheme, where the couplings flow with the renormalization scale ji. The 1-loop RG equations for 
the Higgs-Yukawa model are 

d7i 1 

11— = —^{W^ + mpXy^-AmFy'^). (3.6) 
a}X Ibn^ 

We can set the initial conditions = v) = 3mjj/v^ and = v) = mt/v. If rut is sufficiently 
heavy relative to mh, the Yukawa coupling dominates the RG flow and dX/djl < 0. The renor- 
malized Higgs coupling eventually becomes negative at some pL = E. If the instability occurs at 
very large <p /v, large logarithmic terms ln(0 /v) in U^n might spoil the perturbative expansion. This 
can be reduced using renormalization group improvement to resum the leading large logarithms. 
The dominant terms of U^n at large then become A(jU)(/>^(/x). Hence X{E) =0 indicates that the 
ground state is just about to become unstable. 
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3.3 The constraint effective potential on the lattice 

We can calculate the exact effective potential non-perturbatively, using lattice simulations. 



This was first shown in the pure Higgs theory by Kuti and Shen [|15|]. There is some finite lattice 
spacing a on the lattice which restricts the momenta \p^\ < 71/ a replacing the sharp momentum 
cutoff used in section 2. For a Higgs-Yukawa theory with Nf fermions, the Euclidean lattice 
partition function is 

Z = n / ^/0oW[Det(D[0o])f^exp(-5[(/)o]), = n / #oWexp(-5eff[0o]) 



{D[(j)o])xy = Y^d^,xy+yo(l>o{x)5xy, (3.7) 

where the partial derivatives are replaced by finite lattice differences. If the integrand is positive- 
definite, it can be interpreted as a probability density and importance sampling (i.e. Monte Carlo 
integration) can be used to calculate expectation values, e.g. (^o). non-perturbatively with the exact 
distribution [Det(D)]^'' exp(— 5). All dimensionful quantities are calculated in units of the lattice 
spacing a. There is a phase diagram in the bare-coupling space mQ,Xo,yQ. The Higgs phase and the 
symmetric phase are separated by a second order transition, where the vev, va, and the masses mna 
and ifita, vanish. Since the vev and masses are non-zero in physical units, the transition corresponds 
to the continuum limit c? — > 0. To make the cutoff A = n/a large, the bare couplings must be tuned 
to be close to the transition line. If we calculate via simulations that e.g. av = (a^) 0.05 for some 
choice of bare couplings, we can use v = 246 GeV to convert this into a cutoff A w 15 TeV, as well 
as determine mh and nit in physical units. 



In a finite space-time volume D., we will use the constraint effective potential [[Ij, [1^. For a 
pure scalar field theory, this is 

expi-aUam) = n/ [^--^L^i^)) ^M-Sm. (3.8) 

The delta function enforces the constraint that the scalar field ^ fluctuates around a fixed average 
<I>. The constraint effective potential Uq,{^) has a very physical interpretation. If the constraint is 
not imposed, the probability that the system generates a configuration where the average field takes 
the value is 

P{^) = Uxp{-aUam, Z = j d^'txp{-QUa{^')). (3.9) 

This is in very close analogy to the probability distribution for the magnetization in a spin system. 
The scalar expectation value v = (0) is the value of <I> for which Uq, has an absolute minimum. In a 
finite volume, the constraint effective potential is non-convex and can have multiple local minima 
[jlTll. The standard effective potential ?7eff(<J*) is always convex, even in a finite volume, as the 
Maxwell construction connects the various minima. The two effective potentials are identical in 
the infinite- volume limit, \mi£i-^ooUQ,{(^) = ?7eff(<I>), and the constraint effective potential recovers 
the convexity property. In a finite volume, it is more useful to work with the constraint effective 
potential, where multiple minima can be observed and the transition between the Higgs and sym- 
metric phases is clear. It is also more natural, as the probability distribution P{^) can be directly 
observed in lattice simulations. For the rest of this paper, we drop the subscript Q.. 
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3.4 Hybrid Monte Carlo algorithm and the effective potential 

One way to extract the effective potential from lattice simulations is to generate the ensemble 
of configurations, calculate the average scalar field for each configuration and hence the proba- 
bility distribution P{^). The effective potential is extracted by numerically fitting [/eff(<I>) to P{^) 
using Equation (^!9|). This gives the effective potential for all <I> from one simulation, but with 
limited accuracy. An alternative method is calculate the derivative of the effective potential. For 
the Higgs-Yukawa model with Nf degenerate fermions, the derivative is 

^=m^^ + \x{<^^)<i,-NFy{W)>i>, (w)<j. = (Tr(Z)[0]-i )).!,. (3.10) 

The expectation values (...)<j) mean that, in the lattice simulations, the scalar field fluctuates around 
some fixed average value 0. This method determines the effective potential with greater accuracy 
than fitting the distribution P{^), but the drawback is that a separate lattice simulation has to be 
run for every value of This is the method we use in our investigation of the vacuum instability. 



In this section we use staggered fermions [18, 19], one flavor of which corresponds to four 



fermion flavors in the continuum. With one staggered fermion, the determinant Det(D) is real 
but can be negative due to fluctuations. Then the partition function integrand is not positive- 
definite and Monte Carlo integration cannot be applied. To overcome this problem, we simulate two 
staggered fermions, corresponding to eight continuum flavors, as [Det(D)]^ guarantees a positive- 
definite density. We used staggered fermions only in the very early phase of our simulations. The 
complicated taste structure of staggered fermions with the related rooting issues and the lack of full 
chiral symmetry motivated the switch to chiral overlap fermions which are used now exclusively 
in our Higgs project. Staggered results for the effective potential, which are used here mainly for 
simplicity and pedagogy, have been replaced by simulations with chiral overlap fermions. 

Configurations are generated using the Hybrid Monte Carlo algorithm pO|], where a fictitious 
time t and momenta 7i{x,t) are introduced. New configurations are generated from the equations 
of motion 



= K{x,t) 



7C{x,t) 



dSsff 1 dSss 



where the effective action 5eff is given in Equation The second term in 7c{x,t) is included to 
enforce the constraints 

^Y.^{y,t)=^, Y^niy,t)=0. (3.12) 

y y 

We work with fixed lattice volumes of size 8^* x 16. The scalar field has periodic boundary condi- 
tions, the fermionic field is periodic in the short directions and antiperiodic in the long direction. 
We use the standard leapfrog method to solve the equations of motion, where the step-size At is 
adjusted to achieve acceptance rates well above 90%, and each trajectory length satisfies NfAt > 1. 
For each simulation, we generate at least lO'* configurations and check that correlations between 
the configurations are small. 

The basic quantities of the theory are the bare fields and couplings. A particular choice of 
bare couplings puts us somewhere in the phase diagram and all physical quantities are now fixed. 
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A separate constrained simulation is run for each value of to calculate the effective potential 
derivative. The expectation values we measure on the lattice are bare ones, so the simulations give 
the bare equivalent of Equation ( 3.10| ), namely 



which is converted using the relationship between the bare and renormalized fields, 



d^ 



d<t>() 



(3.13) 



(3.14) 



We measure the wave function renormalization factor in separate unconstrained simulations. 

We want to follow the behavior of ?7eff as we approach the continuum limit, the critical surface 
in the bare-coupling space. We make an arbitrary choice 3^0 = 0.5 and Ao = 0.1. The distance 
from the continuum limit is determined by the remaining bare coupling Mq. We obtained results 
for three choices ml = 0.1,0.25 and 0.29. Typical non-perturbative measurements of the derivative 
dUefi/d^ are shown in Figure ^ All dimensionful quantities are in lattice units, e.g. a • <I>. What do 
we expect to see? In the Higgs phase of the theory, Uen should have a local maximum at the origin 
and a local minimum for some non-zero a • <I>. If the vacuum is stable, the local minimum is in fact 
an absolute one. Let us first look at the results for wiq = 0.1, shown in Figure ||. The simulations 
show that dUcfi/d^ vanishes at the origin and at a • <I> « 2.0; these are the extrema. The derivative is 
negative between these points, so the origin is indeed a local maximum. For a • > 2, the derivative 
is always positive and the local minimum appears to be an absolute one. If the vacuum is unstable, 
dUef[ I d^ should turn negative at large a ■ for which the simulations show no evidence. In these 
units, the lattice cutoff is A = %la and the ratio of cutoff to scalar expectation value is A/v w 1.5. 
This is far from the continuum limit. 




Figure 3: The derivative of the effective potential c/f/gff / for the bare couplings yo = 0.5 , Aq = 0. 1 , = 
0.1, for which the vev is av = 2.035(1). The left side plot is a close-up of the behavior near the origin. 
The circles are the results of the simulations and the curves are given by continuum and lattice renormalized 
perturbation theory. 



We vary the bare mass to get closer to the critical surface and the continuum limit for bare 
masses hiq = 0.25 and 0.29 respectively. The simulations show the same qualitative behavior for 
U^<s- the origin is a local maximum, there is an absolute minimum for some non-zero a<J> and no 



12 



New Higgs physics from the lattice 





ml 


^0 


av = (a^) 


aniH 


anit 


0.5 0.1 


0.1 


0.987(1) 


2.035(1) 


0.521(5) 


0.9977(5) 




0.25 


0.9705(8) 


0.811(1) 


0.297(4) 


0.3906(7) 




0.29 


0.9676(7) 


0.4685(6) 


0.248(3) 


0.2230(3) 



Table 1: The wave function renormalization factor, the renormalized scalar expectation value and the Higgs 
and Top masses, obtained from unconstrained lattice simulations. The bare couplings are those used for the 
lattice measurements of the effective potential f/eff- The estimated errors are in parentheses. 



sign of an instability in the potential. The minimum occurs at £?<!> w 0.81 and 0.47 respectively, for 
which A/v w 3.9 and 6.7, pushing towards the continuum limit. 

The first check of these calculations is to run separate unconstrained simulations with the same 
bare couplings, where L.y0(-^) is allowed to fluctuate freely, and to measure independently v = (0). 
This expectation value should be identical to the value of where U^s has an absolute minimum, 
as determined by the constrained simulations. In the unconstrained simulations, the second term 
for 7i(x,t) in Equation ( 3.11 ) is omitted. The results of the unconstrained simulations are given 
in Table |I} There is indeed perfect agreement between the measurements of (a0) and the loca- 
tion of the Uss minimum obtained from the constrained simulations. The continuum perturbation 
theory calculation of U^n is also shown in Figure ^ We only display the large Np result: not 
surprisingly, for Nf = 8, the Higgs-loop contributions are negligible and can be omitted. We see 
excellent agreement with the non-perturbative simulations for <I> < v, as shown in the left side plot. 
However, the behavior as <I> increases is completely different, as shown in the right side plot. Con- 
tinuum perturbation theory breaks away from the simulation results and predicts that the vacuum 
becomes unstable, with dU^^fi/d^ turning negative. The exact non-perturbative calculation shows 
no indication of this. 

What can we conclude from the comparison? Continuum perturbation theory works well for 
<I> less than and even close to the lattice cutoff A = n/a, as shown by the very good agreement 
with the exact lattice calculations. This is the most that one could have expected. The instability is 
predicted at well above the cutoff, which is completely unphysical and where one cannot expect 
the continuum calculation to apply. The exact effective potential, with the full cutoff dependence, 
is absolutely stable. The standard interpretation of the instability in the continuum ?7eff would be to 
say new physics appears at this energy scale to stabilize the ground state. But the actual cutoff of the 
field theory is far below this scale, especially as we get closer to the continuum limit. The instability 
only appears when the finite cutoff effects are ignored — there is no need for new physics. One 
can ask, is it possible to arrange both the standard ground state and the instability to occur well 
below the regulator cutoff? If so, the instability would be a genuine low-energy prediction. The 
answer is no in the Top-Higgs Yukawa model, if only the standard terms are included in the lattice 
Lagrangian. In this case the only freedom one has is the choice of the bare couplings, and nowhere 
in the coupling-space is a genuine instability seen. If higher dimensional operators are included, 
the Ao > condition perhaps could be relaxed by adding new irrelevant operators, like the ^</>^ 
term, to keep the stability of the cutoff theory intact. This scenario requires further investigation. 

It can be shown in renormalized lattice perturbation theory that the breakdown of continuum 
perturbation theory is due solely to the finite cutoff. A finite cutoff is used in the lattice momentum 
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integrals for the radiative corrections and the counterterms of JJ^s, but otherwise the procedure 
is the same as in the continuum. In Figure |3| we see excellent agreement between simulations, 
and lattice and continuum renormalized perturbation theory for <I>/v < 1. As increases, lattice 
perturbation theory exactly tracks the non-perturbative result, showing a perfectly stable ground 
state. The continuum calculation breaks down, not because of large couplings, but because of the 
neglected finite cutoff. 

4. Wilsonian renormalization group and vacuum instability 

Most of the original work on the consistency of quantum field theory considered only ideal- 
ized theories, supposedly fundamental to describe physics at arbitrarily high energies. Although in 
the previous section on vacuum instability and the related Higgs lower bound problem we found 
a non-removable intrinsic cutoff, the analysis was based on the traditional renormalization proce- 
dure. The Wilsonian viewpoint of the renormahzation group provides a broader and more complete 
perspective on the discussion. 

4.1 Wilson's running Lagrangian 

In the 1970s Wilson developed a new, intuitive way of looking at the renormalization of quan- 
tum field theories based on the flow of effective Lagrangians as generated by renormalization group 
transformations This is based on the realization that physics as we know it seems to be de- 
scribed by effective quantum field theories, which are useful only up to the energy scale Aq where 
new and yet unknown physics is reached. Some smooth intrinsic regularization is introduced (in- 
herited from new UV physics) at Aq which in Euclidean space restricts the length of all four- 
momenta. Physics below the cutoff scale Aq is described by a very general 'bare' Lagrangian 
^(Aq) with an infinite series of local terms, constrained only by symmetries. For any choice of 
the coupling constants in the local terms of the bare Lagrangian, the Euclidean path integral of the 
partition function has to be finite and well defined. The most fundamental constraint on the bare 
Lagrangian is the existence and stability of the functional integral which defines the Euclidean par- 
tition function. If the viewpoint of 'naturalness' is adopted, all the coupling constants of the higher 
dimensional operators are chosen to be of order one in units of Aq. Using Wilson's exact renormal- 
ization group we can consider smoothly lowering the regularization scale to some value say, of 
order the energy scale E far below Aq. To keep physics unchanged, the coupling constants must 
change with the regularization scale. Hence we have a running, or effective Lagrangian (A), 
which flows with A and remains stable at every stage of the procedure in the sense of a convergent 
Euclidean path integral. Since we can use the Lagrangian (Ar) to calculate low energy physics 
at the scale E, it is not the coupling constants at Aq that are important, but those at the scale A^. 
The bare couplings have to be close to a critical surface if mp/j/Ao <^ 1 for the low energy physical 
masses niph of the theory. 

An effective field theory is renormalizable if we can calculate all the S-matrix elements for 
processes with energy scale E, up to small errors which vanish as powers of E /Aq, once we have 
determined a finite number of coupling constants at some renormalization scale Kr ^ E. These 
coupling constants are called relevant; all others are irrelevant. Whatever values we choose for Aq 
(as long as it is large enough) and the irrelevant bare couplings Tj (Aq) (as long as they are natural 
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enough), for a particular choice of the relevant operator set A (A^), the irrelevant operator set T] (A^) 
will be of the order of some power of (A^/Ao). In other words, for any point on the submanifold 
of relevant couplings at A^ there is a flow towards it from a wide variety of initial Lagrangians 
at Ao, all of these being equivalent as far as the values of S-matrix elements for processes with 
energies of order £ A^ are concerned. This more general aspect of renormalizability is the 
concept of universality. An effective quantum theory thus gives us a much more general notion 
of renormalizability than we had in conventional quantum field theory: the regularization need no 
longer be removed, and the irrelevant bare couplings need not be zero. 

It is useful now to adopt the Wilsonian view on the running effective Lagrangian to the Top- 
Higgs Yukawa model we investigated in the previous section. 



4.2 Top-Higgs Yukawa model, vacuum instability, and running Lagrangian 

Adapting the notion of the the running Wilson Lagrangian for the Top-Higgs Yukawa model, 
there are only two marginally irrelevant couplings, X{t) and y{t), in addition to the relevant Higgs 
mass operator. It is important to note that the couplings for increasing t = log(Ao/A) flow from 
bare Xq and 3^0 toward their low energy renormalized values as a function of the energy scale. 
For example, in the large Np limit and for large t values, neglecting the irrelevant couplings, the 



flows are expected to look approximately the same as described by Equation (2.17). The Yukawa 



coupling y{t) will monotonically decrease from its bare value jq towards zero, at the logarithmic 



rate of Equation ( [2.15| ) for large t. The Higgs coupling will start from its bare value A(0) = Ao and 



either it will monotonically decrease, or after some initial rising it will turn around and continue 



to decrease monotonically towards zero, at the logarithmic rate of Equation ( |2.1q ) for large t. In 
the Wilsonian picture, all RG trajectories flow from the general coupling constant space of cutoff 
La grangians ^ {A\))jop-Higgs towards the trajectory specified by ( |2.15[ ) and ( [2. 16[ ) with small but 
calculable corrections from irrelevant operators in the large t limit. 

In the Wilsonian view of the running Lagrangian, the cutoff dependent Higgs mass lower 
bound can be determined in the space of the bare cutoff Lagrangians ^{Ao)Top-Higgs from the 
smallest allowed value of X{Ar) for a fixed Aq/Ar <C 1 ratio where a natural choice for A^ is the 
weak boson mass mz, or the vacuum expectation value v. This calculation is, of course, very hard to 
implement operationally with a large number of bare couplings. The important stability condition is 
the only constraint (with, or without naturalness) on the space of cutoff Lagrangians. For example, 
the choice of Ao < a priori should not be excluded at the cutoff scale Ao, but it requires the 
presence of some positive higher dimensional operator, like X'^^ /A^ -0^, with Aq^^ > 0, to provide 
stability. Whether the Higgs mass lower bound will be necessarily associated with the limit Ao — > 0, 
or the Ao < region also needs to be explored remains an unresolved and interesting question. 

In phenomenological applications an attempt is always made to simplify Wilson's framework 
of dealing with the full space of running Lagrangians. Invoking the Ao/A^ — > limit, only the 
running of the relevant and marginally irrelevant couplings is calculated and the effects of irrele- 
vant operators are ignored. In addition, in the application of RG equations to the vacuum instability 
problem, the simplified equations on A(f) and j(f) are running backward from the mz scale towards 
the cutoff Aq. This interchange of the natural Wilsonian UV IR flow with the IR — > UV inte- 
gration of relevant couplings only is a nontrivial proposition because the Wilsonian RG flow is not 
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known to be reversible, and to set all the irrelevant couplings to zero at the scale Ar = mz would 
require an unknown extension of the space of cutoff Lagrangians ^{Ao)Top-Higgs, if it exists at all. 

In most of the phenomenological RG applications this is not a problem. We believe, however, 
that the RG treatment of the vacuum instability problem requires special care. What corresponds to 
the unstable Ues in Figure || is the running which at some scale to, far below the cutoff scale, 
turns negative as the RG is running backward, from t = logAo/mz towards the cutoff scale t = I. 
It is a signal that higher dimensional operators must play a role to provide a continued stability to 
the theory on all scales. It is unlikely that a positive Aq on the cutoff scale can support this picture, 
forcing the running X{t) to turn positive again and produce an effective potential which will turn 
back positive again after a second minimum which might be lower than the original one where 
the spontaneously broken theory was built (decay of the false vacuum). It is more likely that this 
scenario, if it exists at all, will require the Ao < extension of the space of bare Lagrangians. This 
is an extension which remains largely unexplored and we are just beginning to investigate it. 

4.3 Phenomenology from 2-loop continuum RG 



Vacuum instability was first raised in [ )22| ] and it has since been increasingly refined in ap- 
pUcation to the Standard Model ^ ^ 0, 



33, 34, |, 35, 36]. The 



Running Higgs coupling 



state-of-the-art calculation determines the effective potential to one-loop order, with RG improve- 
ment applied up to two-loop order to the running couplings. 

Results from exhibit the unstable Stan- 
dard Model effective potential for mh = 
52 GeV and m, = 175 GeV, where the in- 
stability appears at = 1 TeV. The lower 
bound shown in Figure |I] is also taken 
from The finite width of the lower 
bound is an estimate of the uncertainty of 
the theoretical calculation, including the 
effect of unconsidered higher-order con- 
tributions. The strict lower bound for the 
Higgs mass can be further refined if one 
allows the ground state to be unstable, 
but demands that the time required to tun- 
nel away from the local minimum at v = 
246 GeV is longer than the lifetime of the 




t=log(mu/GeV) 

Figure 4: The running Higgs coupling is plotted for 
different choices of the Higgs mass from our numerical 
solution of the five coupled 2-loop RG equations for the 
^,y,gi,g2,83 couplings . For input, = 1 75 GeV was used 
with the experimental values of the gi,g2,g3 gauge cou- 
plings. The 1-loop matching of the couplings and the start- 
ing scale of the RG was chosen at mz- 



universe [37, 38 



It is clear that the current experimen- 
tal limits on niH bring the lower bound 
into play. For example, a Higgs boson 
with a mass of 100 GeV should indicate a 
breakdown of the Standard Model around 
50 TeV. However, a Higgs mass in the 
range 160 - 180 GeV apparently allows 
the Standard Model to be valid all the way up to the Planck scale. The occurrence of the vacuum 
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instability mostly relies on the relative magnitudes of and while both renormalized couplings 
can remain small and all three gauge couplings of the SM are included. The perturbative RG 
approach, if cutoff effects can be safely ignored, seems to be on solid footing. However, cutoff 
effects played an important role in Top-Higgs Yukawa models where only the Higgs coupling A 
and Yukawa coupling y drive the dynamics. In this approximation we have shown that vacuum in- 
stability cannot be induced with the SM Higgs potential in the cutoff Lagrangian (the possible role 
of higher dimensional operators to induce vacuum instability remains unclear, as we noted earlier). 
However, in the phenomenological application, all five couplings are running and it is important to 
ask: for the cutoff A at or below the Planck scale Mp, should we expect Top quark induced vacuum 
instability with the SM cutoff Lagrangian without adding new operators? Do we expect a qualita- 
tively different picture when compared to the Top-Higgs Yukawa model? From Figure ^ we find 
that the running A turns negative below the Planck scale for Higgs mass values lower than 135 GeV 
and remains negative when Mp is reached. Further lowering the Higgs mass lowers the scale where 
A turns negative. It remains unclear how these RG flows would be effected by holding Ao > in 
the SM Higgs Lagrangian at some cutoff scale A. How some higher dimensional operators might 
provide a well-defined cutoff theory for the choice Ao < will require further investigation. 



5. Higgs mass lower bound from the lattice 

We would like to outline and implement the first step of a robust strategy to calculate the lower 
Higgs mass bound as a function of the lattice momentum cutoff. The question about breaking 
Euclidean invariance with the lattice cutoff will eventually have to be addressed also. 

5.1 Yukawa couplings of the Top and Bottom quarks 

The third, heaviest generation of quarks consists of the left-handed SU (2) top-bottom doublet 
2l = (^^) and the corresponding right-handed SU (2) singlets Jr, b-^. The complex SU (2) doublet 
Higgs field ^{x) with U (1) hypercharge 7 = 1 is <I> = (^o) where the suffixes -i-,0 characterize the 
electric charge -i-l, of the components. Since 0+ and 0° are complex, we can introduce four real 
components, O = (f^^^^j) and the Higgs potential will have 0(4) symmetry, with broken custodial 
0(3) symmetry, if the Yukawa couplings yt and yh, defined below, are different. The Higgs potential 
in the complex doublet notation has the form, 

1 A 

V{^) = -mW<I>+ — (<I>"^<I>)2. (5.1) 

The Higgs field acquires a vacuum expectation value responsible for the spontaneous electroweak 
symmetry breaking with (^4) = v and the first three components vanishing. The vacuum expec- 
tation value V can be related to the Higgs coupling constant by v = y3/Am// with the relation 



between the Higgs mass mn and m given by Equation ( |2.7[ ). 

Of the four Higgs components three represent Goldstone degrees of freedom, which at finite 
weak gauge coupling become the longitudinal degrees of freedom of the massive weak gauge 
bosons with mass = vg2/2. The fourth component corresponds to the physical Higgs boson 
field. We do not use the Higgs mechanism in the limit of zero weak gauge couplings and keep 
all four Higgs field components where the ^1,^2, ^3 fluctuations represent Goldstone particles with 
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the symmetry breaking in the 04 direction. In the SM Lagrangian all four Higgs components are 
treated on equal footing where ^Yukawa describes the interactions of the SU {2)l doublet Higgs field 
with the quark fields 

■^Yuk^wa =yt -QL^'tR+yb-QL^bR+h.c. (5.2) 

= /T2*J** is the charge conjugate of T2 the second Pauli matrix, yt,yb are the top and bottom 
Yukawa couplings, respectively. When they are equal, the 0(3) custodial symmetry of the Higgs 
potential is preserved after symmetry breaking. For unequal couplings, only the SU {2)l symmetry 
of the Lagrangian is maintained. It is easy to write out the Yukawa couplings in components: 

-S^Yukawa = yt{tL{<l>4 " i<p3)tR + bL{i<t>2 " + (5.3) 

yb{tL{(l>i + ih)bR + bL{ih + ^A)bR] + h.c. 
All masses are proportional to v as they are induced by spontaneous symmetry breaking. 

5.2 One-component Top-Higgs Yukawa model 

We have used lattice simulations to study the Higgs- Yukawa model with a single real scalar 
field coupled to the Top quark with three colors using chiral overlap fermions. This theory has 
only a Higgs particle and no Goldstone bosons, and the Top quark color indices correspond to 
three degenerate fermions. We will not be able to calculate a lower bound directly relevant to 
phenomenology. Our purpose here is to explain in a simpler model how this non-perturbative 
calculation can be applied to a more realistic approximation of the Standard Model. 



The Yukawa interaction Lagrangian in Equation (5.4) has a straightforward chiral lattice im- 
plementation in the overlap formulation where the chiral left-handed and right-handed fermion 
components are precisely defined. The simulation of the full doublet with the heavy Top and much 
lighter b quark would be very difficult on the lattice with two very different mass scales for nit and 
nib after spontaneous symmetry breaking. 

One could choose for a pilot study the degenerate case yt = yb which has a recent lattice 
implementation [ p2| , p3| ]. In this limit, there are three massless Goldstone particles contributing 
to Top-Higgs dynamics. When the weak gauge couplings are turned on, the massless Goldstone 
modes become the longitudinal components of the massive weak gauge bosons via the Higgs- 
Kibble mechanism. The limitation of the four-component model with degenerate quark doublet is 
the artificially enhanced fermion feedback into Higgs dynamics. 

Although the degenerate model of the Top and Bottom quarks is easy to accommodate in our 
Higgs lattice toolbox, we chose the single component Higgs Yukawa model for our pilot study 
with only the Top quark included. When the weak gauge couplings are turned on, one can choose 
unitary gauge to eliminate the three Goldstone components. In this gauge, ignoring the weak gauge 
coupling effects to leading order, one is left with diagonal Top and Bottom quark Yukawa couplings 
where the b quark is decoupled in the yt = limit. This is not a full justification for keeping the 
single Higgs field only, and the price to pay is the absence of feedback from the Goldstone modes 
into Higgs dynamics. Since the primary purpose of the initial phase of our Higgs project is to 
develop a comprehensive Higgs lattice toolbox and test its various uses, the limited one-component 
Higgs field dynamics will provide very useful information. The next logical step will be to restore 
the four components of the Higgs field which requires the b quark, and break the mass degeneracy 
moving toward the yh <^ yt limit. 
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5.3 Phase diagram with chiral overlap fermions 



Lattice Yukawa models with staggered and Wilson fermions were studied before [44, 4^ p^ ]. 
In this work, we adopted the overlap fermion operator to represent the chiral Yukawa coupling 
between the Top quark fermion field and the Higgs field. Although this is the most demanding 
choice for dynamical fermion simulations, staggered and Wilson fermions are not suitable for our 
goals. We discussed some difficulties with staggered fermions in section 3. The difficulties with 
Wilson fermions are worse. It turns out to be impossible to tune to the critical surface of the Top- 
Higgs lattice Yukawa model with Wilson fermions while keeping the Wilson doublers on the cutoff 
scale. This is different from QCD applications of Wilson fermions. 

Our massless overlap Dirac operator is defined as a - D = 1 + 75 sign (H^) with //„ = J^Dw 
where D„, is the usual Wilson-Dirac matrix with a negative mass which for a = 1 has the form 



D„)y, = 35^ - ^£ ((1 + Yn)Un{x-y)6,^y+^ + (1 - Yn)ul,{x)d,,y^^,^ 



(5.4) 



Using the modified 75 = 75(1 — aD) gamma matrix, we define two projection operators, P± 

i = \j/P±,YR,L=P±¥- 

WlYr + ¥r¥l = V^(l - fD) and P{x) = YlYr - 



^{l ±Y5),P± = ^(1±75), and chiral fermion components, Yl^r = VP±,Wr,l = P±V- The scalar 



and pseudoscalar densities are given by S{x) 
\tfRYL = Vr5{^-jD)¥- 
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The Top-Higgs Yukawa model with overlap 
fermions is defined by the Lagrangian 
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1 
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0.14 



%[D + yo-M^-l-D)]w{i , (5.5) 

where the bare overlap fermion field Yo and 
the overlap Dirac operator D were introduced 
earlier. Derivatives are represented by finite 



Figure 5: The vacuum expectation value of the lattice 
field 00 is plotted in lattice spacing units a as a function 
of the hopping parameter for fixed values of Aq = lO^'*, 
yo = 0.35 with 3 colors of the Top quark. The lattice 
size is 12^ x 24 for the plotted data. The complete phase 
diagram can be mapped out by varying Aq and yo to de- 
termine K-£.(Ao Jo)- 



lattice differences in Equation ( |5.5| ) and sum- 
mation over a= 1,2,3 for Top color is under- 
stood. The gauge link matrices are set to the 



unit matrix in Equation (|5.4|). 

The starting point for simulations is the 
phase diagram of the theory in the bare cou- 
pling space of Wq, Ao, and yo. The actual lo- 
cation of the critical surface is determined 



from the condition avo = in a large set of 
non-perturbative lattice simulations. This is shown in Figure |5] where the critical critical hopping 
parameter for a particular choice of bare couplings is calculated. The Higgs part of the lattice 
Lagrangian is parametrized in the simulations as 
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with (^0 = \/2jc0o, and rescaled notation jo = jq^/Ik for the Yukawa coupling. The odd number of 
colors of the single fermion required the application of the Rational Hybrid Monte Carlo (RHMC) 
algorithm for chiral overlap fermion. The first new code we developed was based on This 
is the code which is mostly used in our Top-Higgs-QCD simulations. We also developed a special 
FFT version of the RHMC algorithm which exploited the special structure of the Yukawa coupling 
in the overlap Dirac operator of the Top-Higgs model. In the FFT code, Fourier acceleration is 
used in the evolution of the molecular dynamics trajectories which significantly reduced the auto- 
correlation time between independent configurations. The details of our RHMC algorithms will be 
described elsewhere. 

5.4 Comparison of large Np and Monte-Carlo results 

The algorithm was thoroughly tested in the large A^^ expansion of the model where we sim- 
ulated a sequence of Np fermions, each with 3 colors, which can also be interpreted as the Top 
quark with 3A^f colors. The Np —>■ °° limit of the vacuum expectation value v and the Top mass 
ifit were calculated in rescaled Xq/Nf and yo/^/NF variables for the finite volumes of the sim- 
ulations, for fixed value of wxq. For a particular choice of the rescaled couplings, v and m, are 
plotted in Figure ^ as a function of l/Np. The largest number of fermions was 3Nf = 60 in the 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

1 /Nf 1 /Nf 



Figure 6: The vacuum expectation value v of the scaled Higgs filed is plotted on the left as a function of 
1 /Nf for 3Nf fermion degrees of freedom. The blue dot marks the 1 /Np limit. The right side plot 
shows the Top mass as a function of 1 /Nf with the blue dot marking the calculated 1 /Nf — > limit. The 
lattice size was 12^ X 24 for every simulation point. 

sequence. The solid line indicates the scaled asymptotic value of v and nif. The finite A^^ data 
were numerically fitted with an added 1 /Nf correction term which allows numerical extrapolation 
to the 1 /Nf — > limit with perfect agreement. For example, in the vev test of Figure |6| the fitted 
curve is 1.2562(4) — 0.152(2)/A^f and the large Nf calculation gives 1.2555(7) asymptotically, in 
excellent agreement with the simulations. The sequence of simulations were done with bare param- 
eters yo^Np = 0.7184, Xq - Nf = 10^^, and htq = 0.0637. For the same sequence, the Top quark 
pole mass m, was fitted on the right side of Figure |^ as 0.9727(5) — 0.145{2)/Nf- The inverse 
propagator mass asymptotically is 0.9025 which converts to pole mass m, = 0.9725 at the finite 
lattice spacing a of the simulations by the formula arrit = ln ^+^^ , in perfect agreement between 
simulations and the large Nf prediction. The complete agreement between the analytic large Nf 
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prediction and the Monte-Carlo results provides a very strong cross-check for the correctness of 
our simulation algorithm and the analytic framework. 

5.5 First results on Higgs mass lower bound 

After thorough validation of our algorithm, we turned to a preliminary determination of the 
Higgs mass lower bound in the single component Top-Higgs Yukawa model. The heavy Top quark 
will constrain the lightest possible Higgs for any given cutoff in the single component Top-Higgs 
Yukawa model. The starting point for simulations is the phase diagram of the theory in the bare 
coupling space of wtq, Ao and jq. For every choice of the bare parameter set, the vacuum expectation 
value V and the Higgs and Top masses take some values in lattice cutoff units. Keeping both the 
cutoff and the Top mass fixed in physical vev units, we explore all allowed bare couplings and find 
the lightest Higgs the theory can sustain. Repeating this procedure at various distances from the 
critical surface determines how the Higgs lower bound varies with the cutoff. For the Euclidean 
path integral to exist, we have to require Ao > in the model. We could also consider a more 
general Higgs action where the constraint Ao > is relaxed when positive terms hke are added 
in the higher-dimensional bare coupling constant space of the bare Lagrangian. For now we do not 
include such terms which are part of our ongoing investigations. 

Figure ^ displays our preliminary results which 
are not far from what is expected from the ap- 
plication of the renormalization group. Lattice 
artifacts will require additional interpretation in 
the low momentum cutoff range of the simula- 
tions. 

Adding the QCD gauge coupling 

Our algorithm and simulation code has been ex- 
tended to the Top-Higgs-QCD code of three cou- 
pling constants. The only change is to include 
the SU(3) matrix link variables in the Wilson op- 
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Figure 7: The lowest Higgs mass is plotted as a 
function of the lattice momentum cutoff for three 
different values of the Top mass. All simulation data 
are converted to physical units using v = 246 GeV. 



erator of Equation ( p.4| ) in our construction of 
the chiral overlap operator. The numerical de- 
termination of the phase diagram and the Higgs 
mass lower bound in the extended model with 
X,y,gT, couplings (Top-Higgs-QCD model) is part of our ongoing Higgs project. 



6. Higgs mass upper bound and the heavy Higgs particle 

In this section we will review earlier results on the Higgs mass upper bound from lattice 
calculations and illustrate with the higher derivative (Lee-Wick) extension how a heavy particle 
might be exhibited without contradictions with Electroweak precision data. 

6.1 Higgs sector as an effective field theory 

In the Wilsonian view of section 4, the Standard Model is expected to have some yet unknown 
UV completion above a certain energy threshold Aq. This threshold could be as high as Aq = 
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^Pianck> or as low as Aq = 1 TeV. Below scale Aq the SM is described by the familiar degrees 
of freedom for the known particles, including fermions and gauge bosons, in addition to the four- 
component Higgs field. For illustration, we will choose Aq = Mpianck first in the description of the 
Higgs sector without gauge and Yukawa couplings. Generalization to the full Standard Model does 
not add to the purpose of the discussion here. Lowering the cutoff into the TeV range will be part 
of the discussion. If the Higgs sector is treated as an effective theory, the regulator is chosen for us 
as an intrinsic part of the theory. Euclidean four-momenta are smoothly cut off when their lengths 
exceed some scale Aq. In this way all momentum integrals are made manifestly convergent, and no 
infinities are encountered. The simplest choice is an exponential cutoff function in the propagators. 



Ka{p) =exp 



(6.1) 



which can be built into the Lagrangian (Aq) for non-perturbative calculations. A mass term could 
have been added to in Eq. (6J_) but we simplified the notation for this quaUtative discussion. The 



general 0(4) Higgs Lagrangian at scale Aq = Mpianck is given by 



+ ur ur +T^{rrr + ^d^r adf^r + - , (6.2) 

^Planck '''^Planck '''^Planck 



where summation is implied over a = 1,2,3,4. Only a few higher dimensional operators are in- 
cluded for illustration and the exponential cutoff is implicitly understood in the functional integral 
built on the Lagrangian of Eq. (|6^). 

6.2 Higgs mass upper bound from the lattice 

The highest allowed Higgs mass from the Lagrangian of Eq. (^) was investigated before, 
using lattice cutoff with cgjA^jCg and all other higher dimensional couplings set to zero. Cor- 
rections from the higher dimensional operators are expected to be small, of the order of powers 
of niff/Mpianck unlcss the couplings cg, A^jCs, or any of the other higher dimensional couplings are 
pushed toward asymptotically large values. It is a limit which is considered artificial and far outside 
naturalness bounds. 

Convincing evidence for the Higgs upper bound and its numerical value comes from lattice 
calculations [4^, 50] where the derivatives are replaced by finite lattice differences giving up Eu- 



clidean invariance on the Planck scale. The advantage of the lattice approach is that the full Ao 
range can be scanned from to oo. This is important if the Higgs self-interaction is a marginally 
(logarithmically) irrelevant operator in the triviality scenario. In the limit of infinite cutoff, the 
largest allowed Higgs mass would be driven to zero (triviality of the renormalized Higgs coupling), 
but with the cutoff at the Planck scale we will get a definitive nonvanishing upper bound which is 
saturated at Ao = oo in the lattice approximation. The renormalized Higgs coupling at low energy 
can be defined as the ratio Xr = / where v = 246 GeV is the vev of the Higgs field (the fourth 
component of the 0(4) field), and tur is a renormalized Higgs propagator mass which is related in 
two-loop perturbation theory to the physical Higgs mass by the relation mn = ms[l + '^i^^jz^^r]- 
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Based on non-perturbative lattice studies, we expect that the largest Higgs mass is obtained in the 
A() ^ oo limit. For any choice of Ao in the 0(N) Higgs model we have 

= Mpianck •C(Ao) • (i3iA«)"Prexp( - ^) {l + ^(^i?)}, (6.3) 

with j8i = 5 (A^ + 8) jgp- and k = -\ (3A^ + 14) ^^^^2)2 ■ The relevant choice is N=4 for the Standard 
Model. The non-universal amplitude C(Ao) is determined from matching to lattice calculations in 
the range 2k < Kjmu < 100 pOj leading to the upper bound mu = 145 GeV in the Ao = 0° 
limit, if the cutoff is at the Planck scale. In principle, the lattice cutoff could be replaced by the 
exponential cutoff function of the continuum theory. It would be required to replace the momentum 
square in Eq. ( p] ) by its lattice version and take the inverse lattice spacing much larger than Aq. A 
new amplitude would emerge which could change the numerical value of the upper bound without 
breaking Euclidean invariance at finite cutoff. This is particularly useful when the cutoff is brought 
close to the low energy physical scale. In the discussion of the higher derivative extension of the 
Higgs sector we will show how to insert a heavy continuum cutoff scale in the theory which was 
turned into a practical calculation before [ 11 , |l2 1 . This suggests that the insertion of the exponential 
cutoff scale might be feasible in practical calculations. What remains the most interesting question 
for LHC physics is the lowering the cutoff from the Planck scale into the TeV range. This will be 
illustrated next in the higher derivative extension of the Higgs sector with the scale of new physics 
in the TeV range. 

6.3 Higher derivative (Lee- Wick) Higgs sector 

An interesting extension of the Standard Model Higgs sector was proposed earlier by the 



addition of higher derivative operators using ideas originally discussed by Lee and Wick 1 11, 12, 



51| , [52| ]. Recently a complete Standard Model was constructed on similar principles []13[]. Both 
constructions eliminate fine tuning in the Higgs sector and require ghost particles on the TeV scale 
represented by complex pole pairs in propagators with unusual physical properties. The analysis of 



the heavy Higgs particle from [[111, |12p will be followed in our discussion. 

In the minimal Standai^d Model with SU{2)l x U{1)y gauge symmetry the Higgs sector is 
described by a complex scalar doublet <I> with quartic self-interaction as we discussed in section 5. 



The Higgs potential V{^^<t>), as defined in Equation 5J^, is SU {2)i x U{V)y invariant. It also has 



a global 0(4) w SU {2)l x SU (2)^ symmetry, larger than required by the SU {2)l x U{1)y gauge 
symmetry. Before the weak gauge couplings are switched on, it is convenient to represent the Higgs 
doublet with four real components <p" which transform in the vector representation of 0(4). 

We will include new higher derivative terms in the kinetic part of the 0(4) Higgs Lagrangian, 

1 cos (20) 1 

where summation is implied over a = 1,2,3,4. Also, in this subsection and the next, we use the 
Minkowski metric and a familiar, convenient form of the Higgs potential, V{<^"<^") = — + 
X{<p"<p")^. The higher derivative terms of the Lagrangian in Eq. ( ^^ lead to complex conjugate 
ghost pairs in the spectrum of the Hamilton operator. The complex conjugate pairs of energy 
eigenvalues of the Hamilton operator and the related complex pole pairs in the propagator of the 
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scalar field ghost particles are parametrized by ^ = Me^'®. The absolute value M of the complex 
ghost mass ^ will be set on the TeV scale. The Higgs Lagrangian in Equation ( |6.4| ) describes 
a finite field theory without divergences, or fine tuning. It has a particularly simple form with the 
special choice = 7r/4 of the complex ghost phase, 

■^H = \d^rd^r + 2^ ^d^r nd^r-v{rr) ■ (6.5) 

The — > limit in Eq. ( |6.4| ) requires special attention. In this limit, the ghost particle becomes 
real and to avoid a double real pole in the propagator with problematic behavior, the choice = 
requires to drop the ^-gr Od^<p" Ud^ip" derivative term in the Lagrangian, 

= ^d^rd^r - ^ or - v{rr) , (6-6) 



the starting point of 1 13]. 



6.4 Gauge and Yukawa couplings 



= --:F,yF^^ - -—DF,,DF^^ , (6.7) 



Gauging the Lagrangian ( p.5\ ) remained unpublished before [pSQ. For completeness, we present 
the main results. The construction of the higher derivative U(l) gauge Lagrangian mirrors Eq. (| 
for the special choice = 7r/4, 

1„ _„„ 

4M4 

with U(l) gauge field and F^y = d^By — dyB^. In addition to the massless gauge vector boson, 
the higher derivative term in Eq. ( |6.7[ ) will insert a ghost particle in the spectrum of the Hamiltonian 
with a complex conjugate pole pair parametrized by ^ = Me^'®. For a general complex phase 
an additional term will appear in the Lagrangian, in close analogy with the construction of Eq. (|6!4|). 

The higher derivative Yang-Mills gauge Lagrangian for the SU (2)iy weak gauge field will 
follow a similar construction adding the dimension eight ghost term, 

= -\G",yG"^^ - ^.D^GlyD^G"^^ , (6.8) 
where the notation G"^y = d^jWy — dyW^ + gf^'^W^y/y is used with the covariant derivative D^'^ = 



5"*5n +gf"^'^W!^. Higher derivative Lagrangians, similar to Eq. (|6^, were first introduced by 



Slavnov to regulate Yang-Mills theories []54|]. 

Labeling the components of the complex SU {2)l Higgs-doublet field as <I> = (*o) the gauged 
Higgs sector is described by the Lagrangian ^ = + Jfg + ^Higgs with the Higgs Lagrangian 

^Higg. = {D^^Yd^^+ :^^{D^,D'^D^)\D^D^D^) - V(<D^<D) (6.9) 

where the Higgs potential is V(<I>^<I>) = — i/x^<I>^<I> + A(<I>^<I>)^ and the gauge-covariant derivative 
is D^O = + /|a • + i^B^ <I>. The higher derivative term in the fermion Lagrangian will 
take the form 

^fermion = l^P'V ^ T^Tp^ (6.10) 

Next we will briefly summarize two important features of the higher derivative Higgs sec- 
tor with the ghost mass scale in the TeV range. The RG running of the Higgs coupling freezes 
asymptotically and a much heavier Higgs particle is allowed in extended Higgs dynamics. 
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6.5 Running Higgs coupling in the liigher derivative Higgs sector 



with onset of ghost effects ■ 
running coupling freezes 




Figure 8: Running Higgs coupling in the 
higher Higgs sector. 



This can be illustrated by calculating the scale de- 
pendent one-loop j8 -function within renormalized per- 
turbation theory in the broken phase of the higher deriv- 
ative 0(N) Higgs sector [|5|, In addition to N-1 
massless Goldstone modes, there is a massive Higgs 
excitation and a massive complex conjugate ghost pair 
appears in all N channels, as a consequence of the new 
derivative term in the Lagrangian. On a low energy 
scale /I, when t = log(/x/v) is negative, the j8-function 
is dominated by the Goldstone modes whose one-loop 
contribution is Above the Higgs mass thresh- 

old the massive Higgs loop contribution sets in and 
the j8 -function becomes ^^i^A^(f) which is the familiar 



one-loop form in the minimal mass independent sub- 
traction scheme of the standard 0(N) model. As t in- 
creases, the complex ghost loop becomes increasingly important and well beyond the ghost scale 
M, for t » \og{M /vev), the beta-function will asymptotically vanish. The running coupling con- 
stant X{t) first will grow as t increases, but eventually it will freeze at some asymptotic value 
as shown in Fig. |[ Ghost loops in the higher derivative Higgs model cancel the loops effects 
from the low-energy SM particles in the UV region and this 'anti-screening' effect opens up the 
possibility for such theories to be more strongly interacting than the standard Higgs sector. 



6.6 Scattering amplitudes 

The Higgs particle is defined as the resonance pole in the s-channel Goldstone scattering am- 
plitude. The Goldstone amplitude can be calculated in the higher derivative Higgs sector of the 
0(N) Lagrangian in the large N approximation. In addition, the Higgs particle can be investigated 
directly in lattice simulations of the higher derivative model, just like in the standard Higgs sector. 

In Figure ^ we plotted from [55, 57| the cross section as a function of the ^/s center of mass 



energy in ghost mass units. The location of the complex Goldstone ghost pair in the scattering 
amplitude of the first Riemann sheet is determined by the choice of the phase angle = 7r/4 in 



the Lagrangian of Equation (5.4). The peak in the cross section corresponds to the complex Higgs 
resonance pole on the second sheet of the scattering amplitude. 

Also plotted in Figure ^ is the scattering phase shift as a function of y/s. The phase shift has 
a sharp rise at the Higgs pole; however the cross section and the shape of the phase shift do not 
describe a standard Breit-Wigner shape in the presence of the ghosts and higher derivative Higgs 
dynamics. It is 'unusual' that the phase shift decreases as the energy gets through the real part of 
the ghost mass signaling acausal behavior in the scattering amplitude. It had been argued by Lee 
that this acausal behavior would only occur on microscopic scales, typical of the Compton wave 
length of ghosts, and it will not lead to macroscopic acausal observations. In the large N plots of 
Figure ^, the bare parameters were tuned to tue = 1 TeV for the Higgs mass with the ghost threshold 
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Figure 9: The Goldstone Goldstone scattering cross section and phase shift is plotted against the center of 
mass energy in large-A^ expansion for the PauH-Villars higher derivative 0{N) theory. The input vev value 
is V = 0.07 in M units. The peak corresponds to the Higgs resonance, which is at mn — 0.28 in M units. The 
scattering cross section is completely smooth across the so-called ghost pole locations. 



located at 3.6 TeV. Lattice simulations confirmed similar strongly interacting heavy Higgs physics 
scenarios [l55|, ^]. 



6.7 Heavy Higgs particle and the p -parameter 

We discussed in the introduction that a heavy Higgs particle, beyond the 200 GeV range, is 
not consistent with Electroweak precision data in the perturbative sense. Concerns were raised 
earlier that the heavy Higgs particle of the higher derivative Higgs sector will contribute to the 



Electroweak p-parameter beyond experimentally allowed limits []58|]. Straightforward application 
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Figure 10: Higgs contribution to electroweak vacuum po- 
larization operator. 
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of perturbative loop integrals support this 
concern. However, with new physics on 
the TeV scale (represented by ghost par- 
ticles) the loop integrals are considerably 
different. A crude estimate can be made 
by evaluating the contribution of the vac- 
uum polarization tensors n^jIT^ to the 



'^2<A2 (271)4 {k'^+Ml 



W.tree 



){k^+Ml^J{k^ + ^H{k^)) 



with a sharp momentum cutoff in the TeV range and using the tree level Higgs self-energy operator 
'Lii{k^)- The reduction is quite large in comparison with the 1-loop perturbative formula. Replacing 
the cutoff integral by the Pauli-Villars regulator, which is appropriate for the higher derivative 
theory, we get similar reduction. The effects of the non-perturbative Higgs dynamics represented 
by a complicated r//(^^) operator would have to be determined by non-perturbative simulations. 
If these reduction effects are not sufficient, one might need to add another Higgs doublet to the 
extended Higgs sector in the spirit of recent suggestions p9|]. To exhibit a heavy Higgs particle as 
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abroad resonance, with strong interaction and with acceptable p -parameter, remains an interesting 
challenge for lattice Higgs physics and model building. 
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